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Abstract. Despite the strong evidence that the peaks in 
the spectrum of solar oscillations are asymmetric, most de- 
terminations of mode frequencies have been based on fits 
of symmetric Lorentzian profiles to the Fourier or power 
spectra of oscillation time strings. The systematic errors 
resulting from neglecting the line asymmetry could have 
serious effects on inversions for the solar internal structure 
and rotation. Here we analyse artificial data based on sim- 
ple models of the intrinsic line asymmetry, using GONG 
mode para meters with asymmetries found by one of us 
(Rosenthal |l998b ). By fitting Lorentzians to the result- 
ing spectra, we estimate the likely properties of the errors 
introduced in the frequencies. We discuss whether these 
frequency shifts have a form similar to the near-surface 
layers uncertainties and are therefore suppressed in the 
solar structure inversion. We also estimate directly their 
contribution, if any, in the solar sound-speed and density 
determinations using the SOLA technique. 
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1. Introduction 

There is strong evidence that the observed p rofiles of solar 
oscillation are asymmetric (e.g. Duvall et al. 1993 ). This is 
thought to be a consequence of the localized nature of the 
stochastic driving source, possibly combined with a con- 
tribution due to noise which is correlated with the driv- 
ing (e.g. Gab riel 1994 ; Abrams & Kum ar 1996| ; R oxburg h 
& Vorontsov [l997t Ni gam e t al. |l99^ ; Rosenthal |l998a) 3; 
Nigam & Kosovichev 1998a ) . Yet in most analyses of he- 
lioseismic data the observed Fourier or power spectra have 
been fitted with symmetric Lorentzian profiles. This leads 



to systematic errors in the inferred frequencies, of possi- 
bly quite serious effect on the inversion for, e.g., the solar 
internal structure. 

Just how serious will be the effect on the inversion de- 
pends on the variation of the frequency shift with mode 
frequency {v) and degree (l). Our aim then is first to assess 
how the frequency shift due to fitting asymmetric peaks 
with Lorentzian profiles varies in the / — v plane, and then 
to consider its effect on structural inversion, specifically a 
SOLA inversion for sound speed and density. This contin- 
ues our earlier i nvestig ation of this problem (Christensen- 
Dalsgaard et al. 1998 hereafter Paper I). 

First we need a model of the line asymmetry. We use 
the following simple representation of the asymmetric pro- 
files in oscillation power (see also Rosenthal 1998b): 



Pas(i^) 



Here 
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where At/ is the frequency separation between modes of 
adjacent radial orders. Also, vq is the eigenfrequency of 
the mode, and a.^ is a measure of the asymmetry. Sim- 
ple algebra shows that, provided a\ <^ al, Pas{v) has a 
minimum value of ~ N at x ~ a^. Thus A'^ is a mea- 
sure of the ratio of noise to signal in the power. In the 
limit 04 ^ oo we recover a Lorentzian profile with max- 
imum value unity, and with half width at half maximum 
(measured in terms of x) equal to ai, sitting on a uniform 
background of level N. 

We note that Nigam & Kosovichev (1998a) presented 
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a formally different, but mathematically essentially equiv- 
alent, expression for the asymmetrical line profile. This is 
characterized by the dimensionless asymmetry parameter 
B which, in our notation, is given by i? = —ai/a^. 
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To estimate the systematic error in the frequency de- 
termination we have carried out fits of Lorentzian profiles 
to timestrings of artificial data, assumed stochastically ex- 
cited and with a power envelope (see Section 2 for a more 
precise description) given by Eq. This simulates the 
analysis of the actual data. The fit results in a determi- 
nation of the location x''*''^ of that symmetric Lorentzian 
which best fits the asymmetric power distribution. From 
this, the frequency error S^'^^^v resulting from the assump- 
tion of a Lorentzian profile is obtained as 



Jas) 



-Al 



(3) 



Of course, x'^'^^"^ depends on parameters ai (linewidth), 
a4 (asymmetry parameter) and N (background). In order 
to assess how S^^^^'v varies over the I — v diagram, we 
need to know how linewidth, asymmetry parameter and 
background vary with / and v. We have used GONG data 
to estimate the variation of these quantities. 

Finally, we have used a SOLA inversion technique to 
calculate the apparent differences in sound speed and den- 
sity related to the frequency error J^'^'^'i/. 

2. Predicted shifts using artificial data 

We have considered p-mode power spectra in the vicin- 
ity of a single mode. These were constructed, on the as- 
sumption of stochastic excitation, as normally distributed 
Fourier spectra with zero mean and variance given by 
Eq. (|^). To save time and still have good frequency reso- 
lution, the simulated data frequency resolution was taken 
to be a\ divided by 5. 

Each synthetic power spectrum (P{vi)) was fitted by 
minimizing the negative of the logarithmic likelihood func- 
tion (L): 



^(a)^-ln(X)^^{logKa,..)] + ^} 



(4) 



where w(a, Vi) is the model of the variance of the spectrum, 
determined by the parameters a. The model was taken to 
be a symmetric Lorentzian profile. 
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where A is the amplitude, z/q is the fitted frequency, 7 
is the half width at half maximum and B is the back- 
ground. As in the fits of real data, the usual parameters 
were fitted: the central frequency, and the logarithms of 
the amplitude, linewidth and background noise level. We 
have fitted over intervals of 100 times ai. Fig. |l| shows 
two examples of 'peaks' in such artificial power spectra, 
and the discrepancy between the true line profile and the 
fitted Lorentzian. 

We have u sed d ata obtained by the GONG project 
(e.g. Hill et al. 199S) to estimate the parameters of Eq. (|l]) 




-0.20 



-0.10 



0.00 



0.10 



0.20 



o 




-0.20 



-0.10 



0.00 0.10 0.20 



Fig. 1. Synthetic power spectra on a logarithmic scale 
generated using ai — 0.0045 and N — 0.003, as a function 
of X. The two plots are two different realizations using 
q;4 = 0.1 and 0.3. The heavy continuous lines show the 
power envelopes Pas (cf. Eq. |l|) used to generate the spec- 
tra, and the dashed lines are the fitted Lorentzian profiles. 



used to generate the synthetic spectra corresponding to 
different points in the I — v plane. The data were de- 
termined in the GONG pipeline by peak-bagging (us- 
ing Lorentzian profiles) individual modes in the averaged 
spectra from six 3-month series of observations (GONG 
months 10-27). For each (n, Z) we computed the m- 
averaged frequency, line width (specified in the GONG ta- 
bles as the FWHM), mode power and background power; 
based on these averages, we determined the parameters 
a\ and TV (assuming the limit 04 — > cxd) as well as Ai/ in 
Eqs (|l|) and (||); the resulting a\ and N are illustrated in 
Fig. %■ 
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Fig. 2. Parameters determining the frequency shifts. 
Panel a) shows the dimensionless peak width ai (mea- 
sured in units of the large frequency separation Ai>) . Panel 
b) shows the ratio N between the background noise and 
the peak power. The stars are the extrapolated values for 
/ < 2. For later reference, modes at low frequency with an 
ai larger than the general trend are indicated by crosses. 
The inserts illustrate the dependence of the parameters 
on mode order n, indicated by the symbol types listed in 
panel b). 



Unfortunately, the GONG data contain few modes 
with I = and none with ^ = 1 or 2. To complement 
our mode set and make a more realistic representation of 
inversions to infer solar structure, we have extrapolated 
the parameters ai and N determined through the GONG 
data for I < 7 modes to estimate those for I < 3. The esti- 
mation was made for the data set ob serve d by the BiSON 
group and described by Basu et al. ( |1997| ). 



X 2500 



Fig. 3. Muhiplets for which GONG or BiSON (for / < 2) 
data have used in this investigation. The crosses corre- 
spond to the crosses in Fig. 0. 



We restrict the analysis to modes with frequency small- 
er than 3510 /iHz; the mode set contains modes of degree 
up to 150 (Fig. |). 

The dimensionless asymmetry parameter ai is deter- 
mined by conditions near the upper turning point of the 
mode, including the excitation and possibly the effects of 
correlated noise. Although the details of these processes 
are as yet poorly understood, it is perhaps not unreason- 
able to expect that they depend little on degree, at least 
for the relatively modest degrees considered here; for in 
this case the horizontal scale of the modes far exceeds the 
scales of the relevant convective processes. Thus we might 
expect ai to be a function of frequency alone. We as- 
sume that this is precisely the case and estimate a^ from 
the dimensional results for low-degree modes in Fig. 4 of 
Rosenthal ( 1998b| ), scaling with the value of /S.u for low- 
degree modes: log(a4) = —2.70-1-9.52 x 10~^ where v is 
given in /iHz. We note that this a^ and our estimate of ai 
are essentially consistent with the asymmetry para meter 
B obtained from analysis of MDI (Toutain et al. 1998 ) 
and GOLF (Thiery et al. 1999) data. 

We estimate the systematic frequency shifts by fitting 
symmetric Lorentzians to asymmetric profiles, using the 
parameter values estimated from the GONG data. For 
each mode, 1000 simulations were performed; averages of 
the estimated dimensionless frequency shifts x'^^^'^ /tt are 
shown in Fig. ^. 

To a considerable extent, x'^^^^ seems to be just a func- 
tion of frequency. This simple behaviour becomes more 
understandable if we consider the dependence of the rel- 



4 



Rabello-Soares et al.: Effects of line asymmetries 




Fig. 4. Dimensionless frequency shifts a;*^**^-' /tt (i.e., J^'^'^'i^, 
in units of the large frequency separation Ai/) obtained 
by fitting Lorentzians to the asymmetric profile Pusiv) of 
Eq. (|l|), using modal parameters based on GONG data. 
Each point corresponds to an average of 1000 simulations. 



Fig. 5. Relative frequency error resulting from the as- 
sumption of a Lorentzian profile, multiplied by the nor- 
malized mode inertia Qni . The continuous curve is a fit to 
the points of a polynomial of degree 6. The crosses corre- 
spond to the crosses in Fig. ||. 



evant parameters on / and i>. As already stated, it is as- 
sumed that the dimensionless asymmetry parameter is 
a function of frequency alone. Fig. || shows that the dimen- 
sionless line- width ai and the noise-to-signal ratio N for 
the GONG mode set are predominantly functions of fre- 
quency. (That this is so for ai was indeed to be expected, 
from the physical properties of the damping.) However, 
they have also some small dependence on order n. 

The dimensional frequency shifts S^^^^'iy will have a 
much stronger Z-dependence, dominated by the depen- 
dence of Av on I. However, it was shown in Paper I 
that Av is essentially inversely proportional to mode iner- 
tia (e.g. Christensen-Dalsgaard & Berthomieu 1991); this, 
therefore, is also the dominant ^-dependence of the S^'^^^v. 



3. Effect on Inferred Solar Structure 

To test the effect on a solar structure inversion (e.g. for 
sound speed or density) of the systematic error from fitting 
asymmetric peaks with Lorentzian profiles, we have run an 
inversion that assumes the frequency differences between 



observation and a reference model to arise from differences 
in internal structure, plus a possible contribution from the 
surface layers: 



Vnl 



^JSM^M dr -I- 



(6) 



We use a SOLA inversion which seeks to estimate the 
sound-speed (or density) difference as a fun ction of po- 
sition within the Sun (see Basu et al. 1996 for details). 
Since the whole inversion process is linear, we do not need 
to create a full set of artificial data: rather, we simply 
use 5'^^^'' Vnl / Vnl as our data, and this corresponds to the 
systematic error that would be added to a real inversion 
of the GONG data due to fitting asymmetric peaks with 
Lorentzian profiles. 

Note that our formulation in Eq. (^) assumes that the 
difference between the observed and model frequencies in 
general contains a contribution of the form J^surfii^ni)/Qni 



where J^^, 



which as indicated is a function of fre- 
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Fig. 6. Relative difference in: (a) sound speed and (b) 
density due to systematic errors in the inferred frequency 
resulting from the assumption of a Lorentzian profile 
(Fig. H). The mode set for the inversions in each panel 
is the same; the inversions assumed the errors of the ob- 
served data in the combined GONG and BiSON data set. 
The dotted line is for the average over 21 + 1 estimations of 
(5(as)j^ (see below), whereas the dashed line used the J^'^^V 
calculated with MDI linewidths (see below). 



quency alone, is determined by the near-surface errors in 
the model, and Qni is the mode inertia normalized by 
the inertia of a radial mode of the same frequency (e.g. 
Christensen-Dalsgaard & Berthomieu 1991). The SOLA 
inversion allows us to suppress the contribution from the 
function .7-"surf by imposing a constraint that would com- 
pletely remove such a function if it were a polynomial of 
degree A or smaller. Typically we choose A = 6: alter- 
natively we can simply choose not to impose this extra 
constraint at all. 

If the frequency shifts introduced by mis-fitting 
the asymmetric peaks were also to be of the form 
^sm:f{i^ni)/Qnh they too would therefore be suppressed 



Fig. 7. Relative difference in: (a) sound speed and (b) 
density between Sun and the reference model using the 
observed frequencies calculated by the GONG and BiSON 
projects (Lorentzian fitting). The dotted line is the solu- 
tion has been obtained by subtracting the inferred rela- 
tive differences (circles) in Fig. ^, i.e., it shows the solu- 
tion 'corrected' for the effect of fitting the data with a 
Lorentzian profile. 



in the inversions. In Fig. ^, we have plotted the rela- 
tive frequency error resulting from the assumption of a 
Lorentzian profile, multiplied by the normalized mode in- 
ertia Qni- As Qni is asymptotically related to Ai/, Fig. || 
is very similar to Fig. ^ The scattered points at low fre- 
quency (crosses) are associated with the scatter in ai 
(crosses in Fig. |2|). The linewidth determination at low 
frequency is particularly difficult because of its small mag- 
nitude. We have looked at other data sets, e.g. MDI data, 
and this scatter at low frequency is not present, which 
suggests that it is an artifact of a poor determination (see 
Fig. ||a, below). 

Fig. 1^ shows that the relative frequency error due to fit- 
ting Lorentzians, when scaled by mode inertia, is in fact 
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largely just a function of frequency; thus one might ex- 
pect that much of it would be suppressed by the sup- 
pression of the surface term in the SOLA inversion. The 
actual inferred sound-speed and density differences from 
performing SOLA inversions of S^^^'^i'ni are illustrated in 
Fig. ^. For comparison, we show also in Fig. ^ the in- 
ferred sound-speed and density differences between the 
Sun and Model S of Christensen-Dalsgaard et al. ( |1996D , 
using the m-averaged mode frequencies obtained by the 
GONG project plus the BiSON mode frequencies for / < 3 
(described in Section 2). For consistency with the results 
for the solar data, in the inversions we assumed that the 
errors in the S'^'^^Vni were the same as the errors in these 
observed frequencies. The important conclusion is that the 
change in the inferred sound speed due to S'^'^^'^Vni is small 
compared with the total sound-speed difference between 
the Sun and the model. The same is true of the change in 
inferred density below the convection zone, for r < 0.7 R. 
However, it appears that the inference of density in the 
convection zone is more sensitive to the effects of asym- 
metry, as described here. The magnitude of the effect can 
perhaps be better appreciated in Fig. 0. There the dot- 
ted lines show the result that would be obtained from 
inversion of solar frequencies corrected for the effect of 
asymmetry, i.e., from observed frequencies fitted with an 
asymmetric profile given by Eq. (|l|), on the assumption 
that our representation of asymmetry is correct. The only 
significant modification from the original solar inference is 
for 6p/ pair > 0.7 R: there we see that the inferred density 
difference in the convection zone is no longer nearly con- 
stant, in contrast with our usual experience with density 
differences betw een pairs of solar models (e.g. Christensen- 
Dalsgaard |l99|). 

We note that it actually makes rather little difference 
in this case whether the SOLA inversion explicitly sup- 
presses surface contributions or not: the variation with 
frequency in Fig. |^ has little effect on the inversion even 
when the surface constraint is not applied. In fact, even 
without taking explicitly account of a surface term, the in- 
version provides some suppression of contributions of this 
form: this follows from the fact that the near-surface ef- 
fect of, e.g., 8(? jf? in Eq. (^ is also of a form similar 
to the surface term; thus the localization implicit in the 
inversion leads to a partial elimination of such contribu- 
tions. (Of course, the surface term in frequency differences 
between the Sun and adiabatic frequencies of a model is 
typically more than an order of magnitude larger than the 
differences in Fig. ||, so it is indeed important to apply a 
strict surface constraint when inverting real data.) 

It follows that the results of the inversions, shown in 
Fig. ^, are dominated by the scatter of the points in Fig. || 
around the curve fitted as a function of frequency. Much 
of this scatter could be due to noise in the GONG pa- 
rameter values. This then affects the inversions in much 
the same way as random noise in the frequency data. We 
have carried out inversions of normally distributed ran- 
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Fig. 8. Relative difference in: a) sound speed and b) den- 
sity due to 100 different noise realizations (see text). The 
thick line is the relative difference due to same as 

the circles in Fig. ^. 

dom numbers with variances corresponding to the spread 
in the points in Fig. ^ from the fitted curve (neglecting the 
low- frequency values shown by crosses), to demonstrate 
the effect of the scatter on our inversions. The results are 
shown in Fig. ||, indicating that the sound-speed results 
are essentially consistent with such a random distribution. 
In addition, we have verified that the monotonic increase 
with r in Jp/p in the convection zone results from the scat- 
tered points at low frequency, shown as crosses in Fig. ^. 

4. Discussion 

As described in Section 2, when fitting the observations, 
the mode parameters were estimated for each (ti, ^, m) and 
averaged over the (2Z-I-1) different m to give the values for 
the multiplet (n, X). To reproduce the properties of the ob- 
servations, it would therefore be more realistic to perform 
2Z-I-1 simulations and average the estimated b'^'^^v., instead 
of averaging 1000 realizations as considered so far. This 
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will increase the scatter in the frequency shift, especially 
for low-degree modes. The dotted lines in Fig. ^ show the 
corresponding sound-speed and density differences. The 
monotonic decrease in sound-speed and increase in den- 
sity differences towards the surface are still present, some- 
what enhanced for sound speed as a result of the larger 
scatter in S^'^^^v. However, the variations are still generally 
small compared with the total difference between Sun and 
model (Fig. 0). 

The inferred frequency shifts depend rather sensitively 
on the parameters assumed for the modes, particularly 
their line widths, and hence the use of parameter values 
based on just a single data set might be cause for some con- 
cern. We have already noted that noise in the GONG pa- 
rameter determination might have a substantial influence 
on the inversion results. Another point is that the GONG 
pipeline may systematically overestimate the linewidths. 
Evidence for this assertion is shown in Fig. ||a, comparing 
the GONG linewidths used here with MDI linewidths for 
a period corresponding to GONG months 29-31 obtained 
by the MDI Medium / Program (Schou 1998); evidently, 
the MDI linewidths are systematically smaller than those 
from GONG. It seems probable that the systematic dis- 
crepancy arises from the fact that the MDI parameter 
determination takes the m-leakage into account, whereas 
the GONG determination did not. (As an aside, we also 
note that the scatter at low frequency in the GONG data 
- crosses in Fig. |^a ~ is not present in the MDI data.) 

On the basis of this comparison it would obviously be 
of interest to repeat our study using parameters obtained 
from the MDI pipeline. Unfortunately, there is not at the 
moment a good determination of the background power 
when using the MDI pipeline (the fit is performed only 
in a very narrow window around each peak, so that the 
background far from the peak cannot be determined re- 
liably). However, we have tried using MDI linewidths to 
estimate ai and GONG data to estimate A^: the resulting 
f are shown in Fig. (gray dots) and compared 
with the results we obtained with GONG linewidths (black 
dots). As one would expect, the smaller MDI linewidths 
result in smaller values of S^'^^^v. Thus, the results in the 
present study may be pessimistic in the sense that if the 
GONG linewidths are overestimated, then this will cause 
us also to overestimate the effects of line asymmetry. Be- 
sides, the MDI linewidths do not show a dependence on 
degree (or order) at given frequency, in contrast to GONG 
data (cf. Fig. ||a). Thus the dependence of d^^'^^v on de- 
gree is small; the degree dependence that remains is prob- 
ably introduced mostly by the N parameter obtained from 
GONG. The dashed lines in Fig. ^ show the resulting in- 
ferred sound-speed and density differences. The monotonic 
variation towards the surface is substantially smaller than 
for the full GONG parameters, particularly for density, 
and has the opposite sign; this confirms our impression 
that the inferred structural variations in the solar inte- 
rior (Fig. |6[ circles) are probably largely an artifact of the 
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Fig. 9. a) Normalized linewidth obtained from GONG 
(black dots) and from MDI (gray dots) for a period cor- 
responding to GONG months 29-31 by Schou ( |1998[ ). b) 
Relative frequency error resulting from the assumption of 
a Lorentzian profile multiplied by the normalized mode 
inertia Q„; using linewidths obtained from: GONG (black 
dots) and MDI (gray dots). 



scatter in the GONG parameter determinations, which we 
have then used in our asymmetry model. 

5. Conclusions 

Our principal conclusion is that the frequency shift that 
arises from erroneously fitting the asymmetric line pro- 
files with symmetric Lorentzians (as is commonly done at 
present) is rather benign: it is predominantly of the same 
form as a structural near-surface contribution - viz. of the 
form J-sui:iivni)/Qnh ^ function of frequency divided by 
mode inertia, which is largely suppressed in the inversion, 
and much of the departure from this behaviour is likely 
to be caused by observational scatter in the GONG mode 
parameters which form the basis for our study. Thus, al- 
though there is some residual error introduced into the 
sound-speed and density inversions at depth, its magni- 
tude is generally very small. We therefore find no evidence 
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to suggest that ignoring line asymmetry has compromised 
the hehoseismic structural inversions published to date. 

One might worry that since the values of N and ai ob- 
tained from the GONG tables are themselves the result of 
a fit of a symmetric Lorentzian profile, they are also sub- 
ject to systematic error. The symmetric fit picks out the 
power level far from the peak as being the noise level. This 
systematically over-estimates the true noise level which is 
given by the power minimum or trough which lies close 
to the peak. However, that error depends only on ai and 
a4 and should therefore also be a function only of fre- 
quency. On the other hand, the amplitude and line width 
are only moderately affected by the asymmetry (see Paper 
I). Hence we do not expect this to affect our conclusions 
that x^^^") is predominantly a function of frequency and 
that S^'^^^u has predominantly the same functional form 
as a near-surface contribution. 

A recently published letter by Toutain et al. (1998) 



purports to show a much more significant effect on the 
inferred sound speed in the solar core from neglecting line 
asymmetry when determining mode frequencies. Although 
we have not considered precisely the same mode set that 
they did, we view the result of Toutain et al. with some 
caution. They compared the inversion of mode frequen- 
cies obtained by fitting observational data with symmet- 
ric Lorentzians with the inversion of mode frequencies in 
which the low-degree modes only had been fitted with 
asymmetric profiles. By fitting the low degrees asymmetri- 
cally, and the rest symmetrically, it is quite probable that 
one will artificially introduce an Z-dependent error that 
does not scale like inverse mode mass and which will be 
erroneously interpreted in the inversion as a spatial vari- 
ation of the structure in the solar core. 

We did a similar experiment with our artificial data. 
In Fig. we compare the inversion of mode frequen- 
cies obtained by fitting observational data with symmet- 
ric Lorentzians (triangles) with the inversion of mode fre- 
quencies which for the low-degree modes (/ < 2) only have 
been corrected for the effect of asymmetry, by applying 
our estimates of the frequency shift (circles with error 
bars). Note that the inversion of the symmetrical fits (tri- 
angles) and the inversion of mode frequencies which have 
all been corrected for the effect of asymmetry (squares) 
agree quite well in the solar core, as already shown in 
Fig. I^a. However, using the asymmetrical fits only for / < 2 
modes introduces a change in the solution in the solar core; 
this confirms our concern that such an inconsistent treat- 
ment of asymmetry introduces an artificial Z-dependence 
and hence may affect the inversion in the core. We note 
that the effect obtained here is more modest than that of 



Toutain et al. (1998); the effect we find is quantitatively 
similar to what has been found by S. Turck-Chieze (private 
communication). We have been able to make some qualita- 
tive comparisons with the asymmetry corrections applied 
by Toutain et al., thanks to T. Toutain (private commu- 
nication) . It appears that their asymmetry corrections for 
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Fig. 10. The circles with error bars show the results 
obtained with asymmetrically fitted low-degree {I < 2) 
modes and symmetrically fitted otherwise {I > 2). The tri- 
angles represent the result of the symmetrical (Lorentzian) 
fits to all modes (same as circles in Fig. ^) and the squares 
the asymmetrical fits to all modes (same as dotted line in 
Fig. [T^ ). To be compared with Fig. 5 in Toutain et al. 
(|T99M. 



I = 1 and / = 2 are very similar and, encouragingly, these 
agree with our simulated S^^^^iy quite well; but the I = 
corrections are rather different. Evidently such an effect 
would introduce some degree dependence and therefore 
radial variation into Toutain et al.'s inversion results. 

Although we conclude that ignoring line asymmetry 
has likely not compromised helioseismic structural inver- 
sions to date, we should like to emphasize the importance 
of taking into account the asymmetric profile in the esti- 
mation of mode parameters now and in the future. With 
the SOHO satellite and the GONG network giving us 
longer time series and with a better signal to noise ratio, 
the effect of asymmetry on the mode parameter determi- 
nations is significant. This is particularly true as we shall 
be looking for ever more subtle features in the solar inte- 
rior. Another important reason for fitting an asymmetric 
profile is the possibility of combining different observables. 
For example, for solar oscillations observed in Doppler ve- 
locity and continuum intensity, the asymmetry in their 
power spectra has an opposite sign (Nigam & Kosovichev 



1998b) which can lead to different estimates of the mode 
parameters when fitting a symmetric profile. An inversion 
of datasets with incompatible frequencies due to different 
line asymmetries will be seriously compromised, because 
the combined frequency error will not look like a single 
near-surface term and will in general be interpreted by 
the inversion method as a spatial variation inside the Sun. 
Likewise, it is possible that non-contemporaneous data or 
data from instruments with different levels of background 
noise would likewise contain systematic errors with serious 
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conscqiicnccs for inversions, unless the peakbagging takes 
line asymmetry into account. We finally note that the de- 
termination of the asymmetry, and other parameters of 
the modes, provides crucial information about the exci- 
tation about the solar oscillations (e.g. Rosenthal 1998a; 
Kumar & Basu 1999; Nigam & Kosovichev 1999). 
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